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The crystal structures of the Iigand-binding domain (LBD) of the 
vitamin D receptor complexed to lo^SfOHfeDs and the 20-epi 
analogs, MC1288 and KH1060, show that the protein conformation 
is identical, conferring a general character to the observation first 
made for retinoic acid receptor (RAR) that, for a given LBD, the 
agonist conformation is unique, the ligands adapting to the bind- 
ing pocket. In all complexes, the A- to D-ring moieties of the ligands 
adopt the same conformation and form identical contacts with the 
protein. Differences are observed only for the 17/3-aliphatic chains 
that adapt their conformation to anchor the 25-hydroxyl group to 
His-305 and His-397. The inverted geometry of the C20 methyl 
group induces different paths of the aliphatic chains. The ligands 
exhibit a low-energy conformation for MC1288 and a more 
strained conformation for the two others. KH1060 compensates 
this energy cost by additional contacts. Based on the present data, 
the explanation of the superagonist effect is to be found in higher 
stability and longer half-life of the active complex, thereby exclud- 
ing different conformations of the ligand binding domain. 

20-epi analogs j nuclear receptors | VDR 

Vitamin D, a seco-steroid hormone (steroid with an opened 
B-ring) is active only in the metabolized form la,25- 
dihydroxyvitamin D3 [la,25(OH) 2 D3], The genomic action of 
this hormone, like other steroid hormones, is mediated by a 
nuclear receptor (NR): the vitamin D nuclear receptor (VDR). 
This receptor exhibits the characteristic modular organization of 
five to six domains of NRs but presents a rather short N-terminal 
A/B domain and an insertion in the ligand binding domain. 
VDR modulates the expression of numerous genes in target 
cells; these genes are involved in calcium metabolism and bone 
formation and in control of cell growth and cell differentiation 
(1. 2). VDR acts by binding to specific DNA sequences as a 
heterodimer with RXR and to the basal transcription machinery 
in a ligand-independent (TFIIB; ref. 3) and -dependent manner 
(TFIIA; ref. 4). ' Other direct contacts have been revealed 
between the receptor and TATA box-binding protein (TBP)- 
associated factors (TAFs; ref. 5), components of general factors 
required for accurate and regulated initiation of RNA polymer- 
ase II. These factors are specific transcriptional coactivators in 
mammalian cells. A different class of proteins has been charac- 
terized that stimulate the transcriptional activity of liganded NRs 
in an activation function 2 (AF2)-dependent way. These coac- 
tivators are thought to bridge the NRs to the transcriptional 
apparatus. In particular, VDR is regulated by coactivators 
belonging to the steroid receptor activator (SRC)/pl60 family of 
proteins, which contain several LXXLJL motifs (6). This motif 
has been shown to form an amphipatic a-helical structure that 
can interact with the AF2 region of NRs (7). Another class of 
coactivators [vitamin D receptor-interacting protein (DRIP), 
thyroid hormone receptor-associated protein (TRAP), activa- 
tor-recruited cofactor (ARC); refs. 8-11] has been isolated as 
multiprotein complexes and strongly potentiated transcription 
mediated by VDR/RXR in a ligand-dependent manner on 
DNA templates assembled into chromatin (8). One of their 



components, DRIP205, interacts directly with the iigand-binding 
domain (LBD) in the presence of ligand and anchors the rest of 
the subunits to the receptor. 

Among the several synthetic analogs of vitamin D, the 20-epi 
compounds, which exhibit an inverted stereochemistry at posi- 
tion 20 in the flexible aliphatic chain, have attracted much 
attention. They are potent growth inhibitors and inducers of cell 
differentiation, while showing an affinity similar to vitamin D for 
VDR (11). KH1060 (Fig. 1), a member of this 20-epi family, 
exhibits similar properties, with decreased calcemic side effects. 
These compounds induce VDR-dependent transcription at con- 
centrations at least 100-fold lower than the natural ligand and 
present antiproliferative activity several orders of magnitude 
higher than the natural ligand (11-13). The differences in 
biological activity of la,25(OH)2D3 and the 20-epi molecules in 
general, and KH1060 in particular, are known to be VDR-LBD 
dependent, but are not yet understood. The ability of 20-epi 
analogs complexed to VDR to induce transcription correlates 
with the efficiency of these compounds to interact with coacti- 
vators (12). The different proteolysis pattern of VDR in the 
presence of the natural ligand and the 20-epi compounds has 
been interpreted as reflecting large conformational changes of 
the receptor on binding of the latter class of molecules. It has 
been suggested that the 20-epi compounds induce an alternate 
conformation of the VDR-LBD, rendering the receptor more 
resistant to protease digestion (8, 11, 12, 14-16). To investigate 
the binding mode of the 20-epi analogs to the VDR-LBD, we 
determined the high resolution crystal structures of the VDR in 
complex with MC1288 and KH1060 (Fig. l)and compared them 
with that of the natural ligand. The comparison provides the 
rational basis to analyze the various contacts between the ligands 
and the receptor. These contacts are most likely responsible for 
the different biological responses. 

Materials and Methods 

The human LBD of VDR (118-427 A166-216) was purified 
following the protocol described by Rochel et al (17). Crystal- 
lization of VDR in the presence of each of the individual ligands 
la,25(OH) 2 D 3 (la,25-dihydroxyvitamin D 3 ), MC1288 (20-epi- 
la,25-dihydroxyvitamin D 3 ), and KH1060 (la ; 25-dihydroxy-20- 
epi-22-oxa-24,26,27-trihomovitamin D 3 ; Fig. 1) was achieved in 
the same conditions as for the natural ligand as previously 
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Fig. 1. Chemical structures of 1«,25(OH) 2 D 3 and the 20-epi analogs KH1060 
and MC1288. (A) 1a.25(OH) 2 D 3 : 1a,25-dihydroxyvitamin D 3 . (B) KH1060: 
1a,25-dihydroxy-20-epi-22-oxa-24,26,27-trihomovitamin D 3 . (C) MCI 288: 20- 
epi-1a,25-dihydroxyvitamin D 3 

described, employing the hanging-drop vapor diffusion method. 
Crystals grew in a starting solution containing 10 mg/ml of 
VDR, a 5-fold molar excess of a specific agonist, 50 mM 
Mes-KOH, 5 mM Tris-HCI, 10 mM DTT, and 0.7 M ammonium 
sulfate equilibrated for 4 days at 4°C vs. a reservoir solution 
containing 100 mM Mes-KOH and 1.4 M ammonium sulfate at 
pH 6.0. The homogeneity of the protein and the formation of the 
complexes between the different ligands and the receptor were 
verified by SDS and native PAGE and electrospray ionization 
mass spectrometry. 

The crystals were mounted in fiber loops and flash cooled in 
liquid nitrogen after cryoprotection with a solution containing 
the reservoir solution plus 30% glycerol and 5% polyethylene 
glycol (PEG) 400. In the case of the MC1288-VDR complex, the 
formation of ice rings explains the partial loss of reflections. The 
diffracting quality of the crystals was not affected. The data 
collection was performed at 100 K at the Deutsches Elektronen 
Synchrotron (DESY) synchrotron facilities in Hamburg (BW7b 
beamline) for the MCI 288- VDR complex and at the European 
Synchrotron Radiation Facility (ESRF) in Grenoble (beamline 
BM30) for the la,25(OH) 2 D 3 -VDR and KH1060-VDR com- 
plexes. For each of the characterized complexes, only one crystal 
was sufficient for obtaining a complete data set at a 1.5-A 
resolution. The crystals were isomorphous and belonged to the 
orthorhombic space group Y2\2{1\ with the unit cells parameters 
as specified in Table 1 . Data were integrated and scaled by using 
the denzo-scalepack package (18). The program cns-solve 
was used throughout the structure determination and refinement 



Table 1. Data collection statistics 



Ligand 


KH1060 


MC1288 


1«,25(OH) 2 D 3 


Space group 




P2,2,2i 




Cell dimensions 








a 


44.49 A 


44.94 A 


44.93 A 


b 


51.87 A 


51.42 A 


51.32 A 


c 


131.39 A 


132.52 A 


132.32 A 


Resolution (last shell) 


20-1.52 A 


15-1.4 A 


20-1.52 A 




(1.56-1.52) 


(1.43-1.4) 


(1.56-1.52) 


Completeness (last shell) 


98.1% (97.3) 


73.4%* (81.1) 


96.2% (66.6) 


flsym (last shell) 


4.8% (29) 


4.3% (23.6) 


5% (12.2) 


<//cr(/)> (last shell) 


20.4 (4.6) 


23.4 (3.7) 


20.7 (9.8) 


Unique reflections 


45,925 


42,311 


46,136 


Redundancy 


4.8 


3.6 


3.4 



*See Materials and Methods. 
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calculations (19). All data were included in the refinement (no 
a cutoff). Initial phase estimates were obtained by omitting the 
la,25(OH) 2 D 3 from the structure of the VDR-la,25(OH) 2 D 3 
complex previously solved at 4°C (17). A rigid body refinement 
was used to correctly place the model into the asymmetric units 
of all of the studied complexes. The phase bias of the initial 
model was reduced by torsion angle molecular dynamics simu- 
lated-annealing refinement (20), followed by Cartesian molec- 
ular dynamics simulated annealing (21) after fitting the ligand 
into the electron density. Subsequent refinement cycles alter- 
nated least-squares minimizations and manual model building by 
using the program O (22). Solvent molecules were then placed 
according to unassigned peaks in the Fourier difference maps. 
All of the refined models showed unambiguous chirality for 
the ligands and no Ramachandran plot outliers according to 
PROCHECK (23). For the VDR-KH1060 complex, the weaker 
density of the electron density map at position C26a of the ligand 
clearly indicates a disorder of this methyl group and therefore the 
ability to adopt different conformations. The final models of 
VDR-la,25(OH) 2 D 3 and VDR-MC1288 complexes contain 255 
residues, with no clear densities for the last four C-terminal 
residues, and also for the residue side chain of His-377. For the 
VDR-KH1060 complex, the final model contains 253 residues 
with no clear density for the two first N-terminal residues, the 
last four C-terminal residues, and residues 375-377 in loop 9-10. 
The volumes of the ligand -binding pockets and ligands were 
calculated as previously (17). 

Results 

Structure Determination. All observed crystal forms are isomor- 
phous, with no significant changes in the protein conformation. 
To obtain crystals suitable for diffraction studies, the different 
ligands and the receptor must be present in the crystallization 
solution. The protein by itself failed to crystallize when incu- 
bated in the same conditions. The crystallization conditions (see 
Materials and Methods) are similar for the three complexes and 
do not represent any harsh change from physiological conditions, 
suggesting that the structural atomic model corresponds to the 
most stable dominant structure in solution. The structures were 
solved by molecular replacement (see Tables 1 and 2 and 
Materials and Methods). High resolution ligand-omit maps cal- 
culated from the refined structure of the three different com- 
plexes allowed to fit the ligands unequivocally (Fig. 2 A-C). 

The atomic models, when compared with VDR- 
la,25(OH) 2 D^ complex, show rms deviations on Ca atoms of 
0.08 A and 0.14 A for VDR-MC1288 complex and VDR-KH1060 
complex, respectively. The rms deviations on all atoms are 0.38 
A and 0.51 A for VDR-MC1288 complex and VDR-KH1060 
complex, respectively. Variations concern only some side chains 
located at the surface of the protein. Contrary to the belief that 
the 20-epi analogs are inducing a different agonist conformation, 
the overall conformation and especially the position of helix H12 
are strictly maintained in all three complexes. Furthermore the 
ligand-binding cavity is unique and conserved for the three 
complexes. The rms deviations of all atoms comprising the ligand 
pocket are 0.09 A for VDR-MC1288 complex and 0.12 A for 
the VDR-KH1060 complex. The sizes of the three ligands are 381 
A 3 , 375 A 3 , and 392 A 3 for la,25(OH) 2 D 3 , MC1288 and KH1060, 
whereas the volume of the ligand pocket remains unaltered in the 
three complexes (660 A 3 ) and the ligands occupy only 57% of the 
volume of the pocket for the VDR-la,25(OH) 2 D 3 and VDR- 
MC1288 complexes and 59% for the VDR-KH1060 complex. 
This value can be compared with the retinoids that occupy 66% 
of the ligand-binding pocket of retinoic acid receptor (17). 

Ligand-Protein Interactions. The interactions between the ligands 
and the receptor involve both hydrophobic contacts and elec- 
trostatic interactions. The A, seco-B, and C/D rings form 
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Table 2. Statistics of structure refinement 



Ligand 


KH1060 


M1288 


1a,25(OH)2D 3 


ffcryst 


21.2% (41,247) 


21.4% (38,055) 


19.5% (41,502) 


Kfree 


23.0% (4,678) 


24.8% (4,256) 


21.6% (4,634) 


rms bond length, k 


0.011 


0.011 


0.005 


rms bond angle, degrees 


1.5 


1.6 


1.2 


Non-hydrogen protein atoms 


1,992 


2,022 


2,022 


Non-hydrogen ligand atoms 


33 


30 


30 


Solvent molecules 


214 


225 


271 


Average B factor protein 


23.5 


14.4 


13.9 


Ligand 


18.9 


12.2 


9.7 


Solvent molecules 


35.3 


27.7 


26.5 



identical contacts as previously described for la,25(OH) 2 D3- 
VDR complex. The hydroxyl groups make the same hydrogen 
bonds, 1-OH with Ser-237 and Arg-274, 3-OH with Tyr-143 and 
Ser-278, and the 25-OH with His-305 and His-397. The deletion 
of 1-OH or 25-OH leads to the largest changes, with a significant 
decrease in binding (1/1000), whereas that of the 3-OH has a 
smaller effect (1/20) (1). This fact is consistent with the 
observation that the Ser-278 Ala mutation decreases the 
transactivation by only 20% (24). The minor transactivation 
decrease observed for this mutation can be explained by the 
presence of Tyr-143, which maintains an anchoring link. A 
mutation of another residue close to the position C4 of the 
A-ring, Cys-288-Ala, results in a severe decrease in transcrip- 
tional activity (8, 24). Cys-288 is involved in numerous Van der 



Waals contacts with neighboring residues (Tyr-147, Phe-150, 
Ser-278, and Tyr-295). Its mutation affects Van der Waals 
contact with the ligand but, more importantly, destabilizes the 
protein architecture of the ligand pocket. 

The specific interactions observed in the three ligand-protein 
complexes involve the hydrophobic contacts of the 17/3-aliphatic 
chains summarized in Table 3. When comparing the natural 
ligand and MC1288, the main difference observed is the posi- 
tioning of the methyl group C21, which results in different 
contacts with Val-300, Leu-309, and His-397 (Table 3). In 
MC1288, the C21 moiety is closer to His-397. Other protein- 
ligand contacts differ, but, to a lesser extent, they involve the 
methyl groups at positions 23 (His-305), 24 (His-397), and 27 
(His-397 and Val-418). In these two complexes, the carbon C22 



25-QH 



3-OH 







Fig. 2. Crystal structures of the VDR LBO complexed to 1a,25(OH) 2 D 3 . MCI 288, and KH1060. Experimental electron density omit map contoured at 2.0 SO of 
{A) 1a,25(OH) 2 D 3 , (fi) KH1060, and (O MC1288. (D) Stereoview of KH1060 bound to the VDR LBD. (£) Close-up view of KH1060 in the ligand-binding pocket. 
Secondary structure features are represented in blue (a-helices) and green (0-strands). The ligand is colored in yellow with the oxygen atoms in red. The volume 
of the cavity is represented in gray. 
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Table 3. Protein/ligand aliphatic side chain distances observed in the 1a,25(OH) 2 D 3 , MC1288, 
and KH1060 VDR-LBD complexes 



Residue 


Atom 


1«,25(OH) 2 D 3 


MC1288 


KH1060 


Leu-227(H3) 


CD1 


C26(3.5) 


C26(3.4) 


C26(3.4) 


Ala-231(H3) 


CA 


— 




C27a(3.7) 


Val-234(H3) 


CB 


— 


— 


C27a(3.7) 




CG1 


— 


— 


C27a(3.8) 




CG2 


C24(3.9) 


C24(3.9) 


C24a(4.1). C27a(4.0) 


Val-300(H6) 


CGI 


C2 1(4.0) 


C2 1(5.4) 


C2K4.2). 022(3.8) 


His-305(L6-7) 


NE2 


C23(3.6) 


C23(3.4) 


C24(3.4) 






C26(3.8) 


C26(3.8) 


C24a(4.1). C26(3.8) 




CD2 


C23(3.8) 


C23(3.4) 


C23(3.9) 






C24C4.6) 


C24(4.5) 


C24(3.6) 


Leu-309(H7) 


CD2 


C2 1(3.8) 


C2 1(4.2) 


C21(3.8) 


His-397(H11) 


NE2 . 


C24(3.7). C27(4.2) 


C24(4.1), C27(4.0) 


C24(4.1). C24a(4.3). C27(4.1) 




CD2 


C21(4.5) 


C21(3.9) 


C2 1(4.4) 


Leu-404(H11) 


CD2 


C26(4.3) 


C26(4.4) 


C26(3.8) 


Val-418(H12) 


CG1 


C27(3.9) 


C27(4.2) 


C27(4.1). C27a(3.9) 



A cutoff of 4.0 A has been used. Underlined contacts in the table are included for comparison purposes. No 
distance is reported for the carbon atom at position 26a; this atomic position is not well resolved (see text). 



makes no contact at a distance closer than 4.2 A. In the case of 
KH1060, the methyl group C21 is quite close to that of the 
corresponding atom in la,25(OH) 2 D 3 . Note that the oxygen 
atom at position 22 of KH1060 forms a van der Waals contact 
with Val-300. The major differences observed between KH1060 
and the two other ligands are the tighter and more numerous 
ligand-protein contacts. The methyl groups C26a and C27a, 
specific to KH1060, form additional contacts with H3, loop 6-7, 
Hll, and H12. A weak density is observed for the C26a methyl 
group, suggesting a structural disorder, whereas the C27a methyl 
group is clearly defined. 

Conformation of the Bound Ligands. In the three complexes, the 
ligands adopt an elongated conformation (Fig. 3) similar to that 
described in the VDR-la,25(0H) 2 D 3 structure (17). The typical 
three ring structure of vitamin D and the position of the three 
hydroxyl groups are preserved among all of the analyzed mol- 
ecules. In all complexes, the ligand is rather tightly bound to the 
receptor around the A-, seco-B, C-, and D-rings. In contrast, the 
aliphatic chain is less constrained, thus allowing alternative 



A 



1-0H 



3-OH 




25-OH 



B 



Torsion angles 


la^5(OH>, 


MC12S8 


KH1060* 


C16-C17.C20-C22 


Jl 


-172 


-164 


C17-C20-C2I(O»)-C23 


-157 


178 


149 


C20-C22(O*)-C23.C24 


164 


•70 


1& 


C22(0")-C23-C24.C2S(24a') 


157 


-168 


-143 


C23-C24-CHa-C25» 






-137 


C23(24»)-C24(24a*)-C25.025 


-60 


.55 


ill 



Fig. 3. {A) Superposition of 1a,25(OH) 2 D 3 (yellow), MC1288 (green), and 
KH 1 060 (blue) ligands after superimposed VDR complexes. Oxygen atoms are 
colored in red. (8) Torsion angles of the 1 70-aliphatic chains of 1a,25(OH) 2 D 3 , 
MCI 288, and KH1060 in the VDR complexes. 



conformations of the 17/3-chain for the natural and the 20-epi 
ligands. The distances between the 25-OH group and the l-OH 
and 3-OH are similar for the three complexes. These represent 
the anchoring points that must be maintained to obtain an active 
conformation. 

In these complexes, the A-ring adopts a chair B conformation 
with the 19-methylene "up" and the la-OH and 3/3-OH groups 
in an equatorial and axial orientation, respectively. Any alter- 
native conformation would affect the H-bond interactions. A 
chair A conformation of the A-ring (la-OH and 3j3-OH groups 
in axial and equatorial orientations respectively) would disrupt 
the hydrogen bond network observed between the hydroxyl 
groups and the protein, and would therefore not fit in the VDR 
ligand-binding pocket because of a steric clash with Phe-150. 
Note that, in all three complexes, a water channel reaches the 
surface of the protein, starting at the position C2 of the A-ring, 
suggesting that a synthetic ligand could occupy this additional 
space without altering the conformation of the A- to D-ring 
conformation. 

A tight channel around the seco-B ring is formed in all of the 
characterized structures by Ser-275 and Trp-286 on one side and 
Leu-233 on the other side. This channel accommodates an 
almost trans conformation of the triene system, connecting ring 
A to the C- and D-rings. The C6 — C7 bond exhibits a torsion 
angle of -149° (C5— Co— CI— C8) that deviates significantly 
from the planar geometry toward the a face of the A-ring and 
confers a bent geometry to the molecule. The a face of the C-ring 
is lined by Trp-286 whereas the methyl C18 on the /3 face points 
toward Val-234. In the case of the VDR-MC1288 complex, the 
CD-rings are slightly tilted, revealing an additional degree of 
freedom to orient the 17/3-chain. 

The remaining section of the pocket is large enough to 
accommodate different variants of the long aliphatic chain of 
vitamin D analogs. The chain flexibility allows the different 
ligands to adapt to the pocket and to form the anchoring 
hydrogen bonds. Fig. 3 shows the dihedral transitions for each 
ligand side chain as observed in the complex. The aliphatic side 
chain of la,25(OH) 2 D^ adopts an extended conformation par- 
allel to the C13— C18 bond with the C13— C17— C2(>— C22 
torsion angle close to 90°. The two 20-epi analogs use different 
strategies to fit into the cavity. KH1060 accommodates its longer 
chain by adopting an eclipsed conformation around the 
022 — C23 bond whereas MCI 288 adopts a gauche conforma- 
tion and has to compensate for a chain that is too short to 
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maintain the 25-OH interactions because of the C20 inverted 
stereochemistry. Our preliminary docking experiments of the 
20-epi analogs have correctly positioned the methyl C21 in the 
same cavity as la,25(OH) 2 D 3 , but failed to give the correct 
geometry of the aliphatic chains. The crystallographic structures 
provide the exact information regarding these changes and the 
details of the interactions with the protein. Compared with the 
natural ligand, the aliphatic chains of the 20-epi analogs are 
lining the opposite side of the pocket. The different orientation 
of the chains is reflected by the different values of the torsion 
angles C16— C17— C20— C22 and CI 7— C20— C22(022 for 
KH1060) — C23 (Fig. 3). Further down the chain, whereas 
la,25(OH) 2 D3 adopts an extended conformation, the MC1288 
and KH1060 exhibit gauche and eclipsed conformations, respec- 
tively, to properly orient the 25 hydroxy! group. The energetically 
unfavorable eclipsed conformation observed in KH1060 
(C20— 022— C23— -C24 of KH1060 equal to 16°) is made pos- 
sible by the strong interactions formed by 022 with CI 2 (2.9 A) 
and C18 (3.24 A). A methylene group at this position as in 
MC1301 would be shifted away from the D-ring. A crude 
modeling analysis suggests that the methylene group would move 
away from CI 2, adopting a gauche (+) conformation instead of 
the eclipsed one observed for KH1060. As a result, additional 
contacts with the protein pockets are formed (VaI-300), which 
may contribute to the higher affinity of MC1301, and other 
contributions like solvation cannot be excluded. Similarly, the 
naturally occurring hydroxylation of C24a in KH1060 and C23 
in la,25(OH)2D3 would form steric clashes with His-305, affect- 
ing the hydrogen bond network with 25-OH. The resulting 
conformational changes can thus explain the reported lack of 
transcriptional activity (25). Based on the crystal structures, this 
activity would depend on the capacity of the pocket to accom- 
modate these additional groups. 

Discussion 

The high resolution structures of the VDR complexed to 
lof,25(OH) 2 D3 and 20-epi analogs provide experimental images 
essential for a rational analysis of the interactions of these analogs 
with VDR and suggest how the 20-epi analogs enhance the stability 
of the receptor and stimulate the interaction with coactivator 
complexes like DRIP. Our data reveal that, contrary to what was 
previously thought, the interaction between the 20-epi agonists and 
the VDR does not result in significant variations from the LBD 
conformational change induced by the natural ligand. The protein 
structures are identical to that observed, in the complex with the 
natural ligand. A similar observation has already been made for 
retinoic acid receptor, where numerous complexes with the two 
natural Iigands (all-trans- and 9-ds-retinoic acid) and synthetic 
analogs revealed that the ligand-binding pocket was unchanged and 
that it is the ligand that adapts (26-29). Such observations should 
be valid not only for retinoic acid receptor and VDR but most likely 
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